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Abstract:  Analytical  solutions  are  derived  for  a  one-dimensional  model  of  the  bulk  temperature  response  of  open-channel  flow  with 
unsteady  and  nonuniform  heating  at  the  upstream  end,  the  water  surface,  and  the  riverbed.  The  solutions  are  explicit  formulas  comprised 
of  transient  terms,  which  play  dominant  roles  in  the  upstream  region,  and  equilibrium  terms,  which  determine  the  temperature  far 
downstream.  The  applicability  of  the  solutions  to  practical  problems  is  illustrated  for  two  cases:  (1)  a  stream  bounded  at  its  upstream  end 
by  a  dam  and  with  a  midreach  inflow,  and  (2)  Boulder  Creek,  Colo.,  which  is  impacted  by  effluent  released  from  a  wastewater  treatment 
plant.  The  model  prediction  is  in  reasonable  agreement  with  gauged  data. 
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Introduction 

Although  accurate  and  computationally  affordable  numerical 
simulations  are  becoming  the  dominant  approach  for  predicting 
mass  and  heat  transfer  in  various  flows,  analytical  approaches 
continue  to  be  useful  because  of  their  irreplaceable  roles.  This 
technical  note  presents  a  model  for  the  bulk  temperature  of  water 
bodies  such  as  rivers,  lakes,  and  reservoirs  that  have  time- 
dependent  thermal  forcing  at  their  upstream  ends  and  water  sur¬ 
face,  side  bank,  and  bed  (hereinafter  referred  to  as  lateral 
boundaries).  The  model  is  based  on  an  earlier  model  described  by 
Edinger  et  al.  (1974)  and  Jobson  and  Schoellhamer  (1987),  and  it 
has  been  applied  in  natural  stream  predictions  (Kim  and  Chapra 
1997;  Boyd  and  Kasper  2003),  variable  reservoir  releases  (Carron 
and  Rajaram  2001),  and  river  temperature  control  (Gu  et  al 
1999).  A  comprehensive  description,  including  supporting  param¬ 
eters  for  phenomena  such  as  air-water  heat  exchange,  can  be 
found  in  Edinger  ct  al.  (1974),  Jobson  and  Schoellhamer  (1987), 
and  Martin  and  McCutcheon  (1999).  Whereas  the  governing 
equations  of  the  model  arc  typically  solved  numerically,  analyti¬ 
cal  solutions  are  sought  in  this  work. 

This  study  assumes  that  the  flow  is  fully  mixed  over  its  cross 
section  with  negligible  effects  of  diffusion  and  that  a  one¬ 
dimensional  model  is  appropriate.  The  temperature  response  is 
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physically  simulated  as  a  scries  of  kinematic  waves  subject  to 
thermal  upstream  and  lateral  boundary  conditions,  which  are  ei¬ 
ther  periodic  or  arbitrary  functions  of  time  and  space.  The  result¬ 
ing  governing  equation  is  a  first-order  linear  differential  equation 
that  can  be  analytically  solved.  The  solutions  extend  the  equilib¬ 
rium  solution  of  Edinger  et  al.  (1974),  which  is  valid  far  from  an 
upstream  region.  Analytical  solutions  have  been  derived  for  simi¬ 
lar  equations  (Chatwin  1973;  van  Genuchten  and  Alves  1982; 
Shukla  2002;  Weigand  2004).  These  solutions  consider  longitudi¬ 
nal  diffusion  and  dispersion  of  flows  with  initial  conditions  and 
unsteady  upstream  boundary  conditions  but  ignore  time- 
dependent  heating  and  cooling  along  the  channels. 


Conceptual  Model  and  Governing  Equations 

The  one-dimensional  governing  equation  for  the  bulk  temperature 
variation  neglecting  diffusion  in  a  semiinfinite  stream  can  be  ex¬ 
pressed  as  (Edinger  et  al.  1974;  Jobson  and  Schoellhamer  1987) 

ST  dT 

—  +  u—  =  K(TM-n  0<5<oo  (1) 

dt  ds 

T  =  Tup,  s  =  0  (2) 

where  /=time,  s=distancc  downstream;  T=  temperature,  Tup 
= temperature  at  the  upstream  end;  =  lateral  boundary  tempera¬ 
ture;  u = flow  velocity;  and  Af=constant  bulk  heat  transfer  coeffi¬ 
cient.  In  order  to  complete  the  general  problem  stated  by  Eqs.  (1) 
and  (2),  an  initial  condition  for  temperature  is  needed.  This  study 
considers  only  equilibrium  solutions  independent  of  the  initial 
temperature  condition,  and  thus  no  initial  condition  is  used. 

Case  1:  Constant  Velocity  Flow  with  Sinusoidal 
Time-Dependent  Thermal  Discharge  and  Boundary 
Heating 

Case  1  has  the  following  flow  and  temperature  conditions: 

u=U  (3) 

JOURNAL  OF  HYDRAULIC  ENGINEERING  ©  ASCE  /  APRIL  2009  /  327 


Downloaded  IS  M  ir  2009  to  134  ?4.20  IS.  Rdlistr  button  subject  to  ASCE  I'censc  01  copyright,  see  http  ''puhs.asce  org/copyi  tyht 


Tup=Tl)+Ti)  sin(o)()/  +  a) 


(4) 

Tm  =  T\ +  TX  sin  wit  (5) 

where  U=  constant  velocity;  T(h  70,  a>0,  and  a  =  average,  ampli¬ 
tude,  frequency,  and  phase  of  the  upstream  temperature  fluctua¬ 
tion,  respectively;  and  Tu  Tx ,  and  a>,  =average,  amplitude,  and 
frequency  of  temperature  at  the  lateral  boundaries,  respectively. 
f0,  r(),  o)0,  a,  Tx,  r„  and  w,  =  constants.  Eq.  (4)  and  (5)  are  fre¬ 
quently  used  as  approximations  of  measurements  (Velz  1970; 
Edinger  et  al.  1974).  For  example,  Eq.  (4)  may  model  the  tem¬ 
perature  modulation  of  water  released  from  detention  ponds  and 
reservoirs,  whereas  Eq.  (5)  may  represent  diurnal  heating  at  the 
water  surface.  Case  I  is  a  direct  extension  of  the  problem  dis¬ 
cussed  by  Edinger  et  al  (1974),  in  which  Eqs.  (I),  (3),  and  (5) 
were  used  to  describe  the  thermal  response  of  natural  water  bod¬ 
ies. 

Case  2:  Flow  with  Spatially  and  Time-Dependent 
Boundary  Heating 

Case  2  is  an  extension  of  Case  I  but  with  a  spatially  variable 
lateral  thermal  boundary  condition  specified  as 


Tm  =  T,  +  T[s  +  (Tj  +  7ji)sin  oij/  (6) 

Eq.  (6)  can  represent  temporal  and  spatial  changes  in  atmospheric 
temperature,  riverbed  temperature,  groundwater  discharge,  and 
snowmelt  at  higher  elevations  (Hanrahan  2007;  Westhoff  et  al 
2007). 

Case  3:  Flow  with  Spatially  Variable  Velocity  and 
Time-Dependent  Upstream  Temperature  and  Boundary 
Heating 

Case  3  incorporates  a  spatially  variable  flow  velocity  and  time- 
dependent  temperature  at  the  upstream  and  lateral  boundaries 


u  =  V(s) 

(7) 

rup=ww 

(8) 

+  2  T,  sin(ay  +  pf) 

(9) 

i-\ 


The  right-hand  side  (RHS)  of  Eq.  (9)  defines  an  arbitrary  function 
of  time  if  the  Fourier  trigonometric  series  converges.  Therefore, 
the  flow  velocity  is  an  arbitrary  function  of  distance,  and  the 
inflow  and  lateral  boundary  temperatures  are  arbitrary  functions 
of  time. 


s/Lt 

Fig.  1.  Base  solution  of  Eq.  (10).  U= 0.1  ms'1,  T|=50°C,  and  K 
=0.04  s*1.  Thj/f^l  +0.4  sin(0.3f),  and  TU[)/f{  =  0.4+0.2  sin(2r 
+  1.5). 


and  the  solution  of  Case  I ,  termed  the  base  solution,  can  be  de¬ 
rived  as 


T=7\  + 


AT, 


\AT2  +  u)^ 


sin(o),r  -  0,) 


ro-T.  +  To  sin 
AT. 


wol'--J  +  a 


7  ■ 

V/f2  +  0), 


“'I'T/  1 


n 

s] 

(J 

exp 

(10) 


where  0i=tan“' (<*>!/ K)  and  LT=Uf  K.  The  base  solution  consists 
of  two  parts:  (1)  the  first  two  terms  on  the  right-hand  side  of  Eq. 
(10),  which  reflect  the  influence  of  the  lateral  boundary  condition, 
are  the  equilibrium  solution  far  downstream  as  discussed  by 
Edinger  et  al.  (1974);  and  (2)  the  exponential  term  on  the  RHS, 
which  results  from  both  the  upstream  and  lateral  boundary  con¬ 
ditions,  represents  a  transient  stage  decaying  with  downstream 
distance.  Lr  is  a  characteristic  transient  length  measuring  the  im¬ 
portance  of  upstream  conditions  The  transient  stage  ends  at  about 
s=3LT  at  which  the  transient  component  will  reduce  to  5%  of  its 
value  at  upstream. 

The  temperature  solution  of  Case  I  oscillates  with  decreasing 
amplitude  until  the  equilibrium  region  is  reached  far  downstream 
(Fig.  I).  The  solution  envelope  is  also  damped  downstream,  but 
the  damping  is  not  monotonic.  This  behavior  can  be  explained  as 
follows.  Let  o)0 =(o,,  Eq  (10)  becomes 


Solutions  and  Discussion 
Case  1 

Eqs.  ( 1 )— (5)  can  be  transferred  into  an  initial  value  problem  for  a 
linear,  first-order,  nonhomogeneous  ordinary  differential  equation. 


T(t,s)  =  T,  +  (T0  -  T,)exp]  — 


+  A  sin(o),/-01  +  k)  (II) 


where 
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A  = 


B  expi  - 


AT,  \2  2BKT | 

K1  Ar  tA2  I  ^  JtP2  Ar  iA2 


I  + 


J 

L, 


+  0I+-Y 


exP|  “ 


j 

z 


(12) 


B  = 


K2  +  i rif 


2A'7~,)7~| 

\KT+Z2 


cos(a  -  0|) 


oscillate  with  distance,  and  the  temperature  envelope  is  defined 
(13)  by  parallel  lines. 


_  Case  2 

T0  sin  a  -  (AT )/  vK2  +  <of)sin  6, 

*Y-tan  The  solution  of  Case  2  is  derived  as 

T{)  cos  a  -  (ATj/VA'2  +  a)|)cos  0, 


,  fl  exp{- j/Lr}sin(- j/L£  + 0!  +  y) 

k  =  tan  - -  - - - 

KTXI\K2  +  o)j  +  B  exp{-  j/Lr}cos(-  s/LE+  Oj  +  y) 

05) 

where  A  =  amplitude  of  the  solution  envelope  and  LE=Uf a),, 
being  a  characteristic  envelope  length.  Eq.  (12)  determines  the 
envelope  wavelength  as  2tt Le.  The  envelope  in  Fig.  I  has  a  spa¬ 
tial  wavelength  of  27TL£/Lr=0.84. 

An  interesting  situation  occurs  in  applications  such  as  river 
temperature  control  when  the  fluctuating  component  of  the  up¬ 
stream  temperature,  Eq.  (4),  is  equal  to  the  equilibrium  tempera- 
ture  far  downstream  as  given  in  Eq.  (10);  T0=KT]I \K2  +  bi],  a)() 
=0)!,  and  a=-0(.  For  example,  using  the  parameters  of  Fig.  I,  the 
resulting  solution  reveals  that  the  envelope  appears  as  parallel 
rather  than  wavy  lines  (Fig.  2)  The  solution  (10)  now  becomes 


T(t,s)  =  T]  +  T\s- 


KUT ; 


,  UT\  K(T\  +  T\s)  . 


K 


sK2 


sin(o),r  -  0,) 


+  a). 


K2  +  o>2 


sin(a> !/  -  0j) 


To“T)+  k  +T()sin,ja)0| 


AT, 


(-5)« 


\K2  +  <o? 

KUT\ 


sin)  | -0 


K2  +  ox) 


sin)  o>|( —  J  - 0J 


exp)  - 


(17) 


—  A  /  |  —  —  I  At 

T(l,s)  =  T |  +  ,  =  sin(a)|/-0|)  +  [r(l-7|]expj  -  —  f 

v'/f2  +  0)?  {  Lt  J 

(16) 

The  coefficient  before  the  exponential  term,  the  third  term  on  the 
RHS,  is  a  constant  and  thus  the  extremes  of  the  solution  will  not 


where  0j  =  tan"l(2a>|A7(A'2-<i>2)).  Eq.  (17)  with  the  flow  velocity 
and  upstream  condition  for  temperature  given  in  Fig.  I,  but  with 
a  temporally  and  spatially  variable  lateral  boundary  heating,  is 
plotted  in  Fig.  3. 


s/Lt 


Fig.  2.  Base  solution  of  Eq.  (10)  with  the  upstream  condition  deter¬ 
mined  by  the  equilibrium  temperature  far  downstream.  U 
=0.11X15-',  f,=50°C,  and  K=().M  s-'.  7^/7,  =0.4+0.0326 

sin(0.3/- 1 .438),  TMIT ,  =  I  +04  sin(0.3/). 


s/Lt 


Fig.  3.  Solution  of  Eq.  (17).  U= 0.1  ms'1,  7',=50°C,  and  K 
=0.04  s'1.  7'bd/7’|  =  1  +0.06.? +(0.4 +0.06.v)sin( 0.3/),  and  rup/f,=  0.4 
+0.2  sin(2r  + 1 .5). 
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Case  3 


The  solution  for  Case  3  is  obtained  as 


T  =  T,  +  2  -  ,  '  ,  sin(uy  +  (3/  -  0;) 
,=i  \  AT2  +  tof 


where  0,=tan  '(a ij/K).  This  solution  will  reduce  to  the  base  solu¬ 
tion  under  boundary  conditions  given  by  Eqs.  (3)-(5).  The  behav¬ 
ior  of  Eq.  (18)  is  illustrated  in  Fig.  4  with  three  examples  that  use 
the  parameters  given  in  Fig.  2,  but  with  different  flow  velocity, 
upstream  temperature,  and  lateral  boundary  condition. 


Applications 

Open-Channel  Flow  with  Time-Dependent  Discharge 
and  Surface  Heating 

Eq.  (17)  is  applied  to  a  river  section  bounded  at  its  upstream  end 
by  a  dam  with  variable  discharge  water  temperature.  The  river 
also  has  diurnal  heating,  which  is  approximated  by  a  lateral 
boundary  condition  that  is  not  only  a  sinusoidal  function  of  time, 
but  also  spatially  variable  because  of  a  change  in  air  temperature 
along  its  length.  The  temperature  variation  at  the  dam  is  specified 
as  having  the  same  frequency  as  that  of  the  surface  heating  but 
with  a  phase  delay.  The  frequency  is  1/24  h  *.  The  river  just 
downstream  from  the  dam  has  a  mean  flow  velocity  of  1  ms'U 
water  depth  of  5  m,  a  mean  temperature  of  15°C,  and  a  daily 
fluctuation  of  3°C.  The  air  has  an  average  temperature  of  20 °C  at 
the  dam  and  25 °C  at  100  km  downstream  with  a  linear  increase, 
an  amplitude  of  15°C,  and  a  phase  lag  of -1.5  rad.  According  to 
Edinger  et  al.  (1974),  K=K/(pCph ),  in  which  p  =  water  density 

(1,000  kg  m-3);  tf=air-water  exchange  coefficient  related  to  fac¬ 
tors  including  wind  speed  (100  W  m"2  °C_1);  Cp- heat  capacity 
of  water  (4,186  J  kg"10^1);  and  h  =  water  depth  (5  m).  Based  on 
these  parameters,  K=  4.8Xl0~6s"1.  The  analytical  solution  is 
plotted  in  Fig.  5,  with  the  upstream  and  lateral  boundary  condi¬ 
tions  given  in  the  figure  caption.  The  envelope  wavelength  2tt Le 
is  86  km.  The  transient  length  3 LT  is  625  km,  suggesting  that  the 
river  section  would  not  approximate  the  equilibrium  within  the 
distance  in  the  figure. 

If  a  multiple  branch  flow  is  considered  with  a  midreach  inflow 
at  5=50  km  that  raises  the  temperature  by  5°C  and  reduces  the 
average  flow  speed  to  0.5  m  s_l,  then  the  solution  for  the  reach 
downstream  of  s  =  50  km  section  is  simply  obtained  from  Eq  (17) 
by  replacing  its  5  with  [5-50]  and  its  upstream  temperature  terms 
(the  terms  with  70)  with  [//(/)  +  5].  Here,  H(t)  =  T(t,5 0),  which  is 
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(a) 


(b)  s/Lr 


Fig.  4.  Solution  of  Eq.  (18).  f/=0.l  ms'1,  T'i  =  50°C,  and  K 
=0.04  s'1:  (a)  u{s)fU=  1  +0.3  sin(3ir$/5),  Tup/ 7\  =04  +  0,2  sin(2/ 
+  1.5),  and  T^!T{=  \  +0.4  sin(0.3/);  (b)  u(s)IU=  1,  7^=04+0.04(1 
+  sin(2/))2,  and  7W=  I +0.4 sin(0.3r);  and  (c)  u(s)IU=  I,  7^/7, 
=04+0.2 sin(2/+  1.5),  and  7^/?!  =  1  +0.4 sin(0.3/)+0,2 sin(0.1r-  1 ). 


the  solution  at  5=50  km  in  case  of  no  inflow  and  directly  deter¬ 
mined  by  Eq.  (17)  The  obtained  solution  is  also  shown  in  Fig.  5, 
from  which  it  is  seen  that  the  solution  is  significantly  changed 
downstream  of  the  midreach  inflow. 
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Fig.  5.  Temperature  response  of  open-channel  flow  with  a  midreach 
inflow,  [/:1ms-1,  rup=15+3sin(7.27X10-s/-l.5)°C,  7^=20 
+  10_\s+ 15  sin(7.27X  IO_5f)°C,  f|=10-5oCm'1,  and  AT=4.8 

X  I0"h  s-1.  /  and  ,s  are  in  seconds  and  meters,  respectively.  An  inflow 
is  located  at  $=50  km,  where  it  raises  temperature  of  the  flow  by  5°C 
and  reduces  the  velocity  to  (7=0.5  m  s-1  for  $5=50  km. 


Creek  Flow  with  Effluent  from  a  Wastewater  Treatment 
Plant 

The  analytical  solution  of  Eq.  (18)  is  applied  to  Boulder  Creek, 
Colo,  which  is  impacted  by  effluent  discharge  from  a  wastewater 
treatment  plant.  A  complete  description  of  the  flow  and  a  compre¬ 
hensive  24  h  survey  can  be  found  in  Windell  et  al.  (1988).  The 
creek  is  13.7  km  long  and  0.2-0.6  m  deep,  the  velocity  ranges 
from  0.12  to  0.4  ms"1.  The  air  temperature,  wind  speed,  solar 
radiation,  and  other  conditions  change  diumally.  The  measure¬ 
ments  of  creek  temperature  were  made  at  four  stations  located 
0.6,  5.0,  9.0,  and  13.7  km  downstream  from  the  plant. 

From  measurement  data  (Windell  et  al.  1988),  w  =  0.27,  3.10, 
0. 13,  and  0.13  ms-1  at  station  0.6,  5.0,  9.0,  and  13.7  km,  respec¬ 
tively,  TM=  15+ 10  sin(7.27  X  10“S/)°C  (t  is  in  seconds);  K 
=  35  W  m~2oC_l;  and  h-OA  m  Here  7^  is  obtained  by  approxi¬ 
mating  the  effects  of  air  temperature  and  solar  radiation.  The 
velocity  is  assumed  to  vary  linearly  between  adjacent  stations, 
and  K  has  been  determined  in  the  same  way  as  in  the  previous 
case.  Letting  Station  0.6  km  be  the  upstream  end  and  approximat¬ 
ing  the  observed  temperature  with  Tup=  17  +  3  sin(7  27  X  10_sr 
-  1.5)°C  [Fig.  6(a)],  Eq.  (18)  is  used  to  predict  the  temperature  at 
the  other  stations,  which  yields  Figs.  6(b-d).  The  analytical  solu¬ 
tion  reasonably  estimates  the  temperature  variations  at  the  other 
stations.  The  largest  difference  between  the  prediction  and  the 
measurement  occurs  at  the  station  furthest  downstream  [Fig. 
6(d)].  Although  the  temperature  predicted  by  the  analytical  solu¬ 
tion  is  not  as  accurate  as  that  from  a  numerical  model  by  (Kim 
and  Chapra  1997),  the  errors  of  the  two  approaches  are  of  the 
same  magnitude.  The  results  arc  promising  given  that  the  flow 
and  temperature  are  affected  by  many  factors  that  the  analytical 
solution  excludes,  such  as  the  material  comprising  the  creek  bed, 
the  thermal  interaction  between  water  and  sediment,  and  the 
three-dimensional  character  of  the  flow.  To  include  these  factors, 
a  numerical  approach  has  to  be  adopted  (Jobson  and  Schoell- 
hamer  1987;  Kim  and  Chapra  1997;  Tang  et  al.  2008). 


(d)  i  (hr) 


Fig.  6.  Temperature  response  of  a  creek  flow  with  effluent 
from  a  wastewater  treatment  plant.  Define  error=[(l  //t)2JL| 
[( T'-Tm)fTm]2]v 2,  where  n  =  total  number  of  measurements  during 
the  24  h,  and  Tm  and  T  are,  respectively,  the  /th  measurement  data 
and  the  corresponding  predicted  temperature,  (a)  Station  located  at 
0.6  km,  error=  0.047.  The  solid  line  is  the  upstream  condition  of  the 
analytical  model  (b)  Station  located  at  5.0  km,  error= 0.098.  (c)  Sta¬ 
tion  located  at  9.0  km,  error=0.056.  (d)  Station  located  at  13  km, 
error=0. 13. 


JOURNAL  OF  HYDRAULIC  ENGINEERING  ©  ASCE  /  APRIL  2009  /  331 


Downloaded  If  Mar  2009  lo  134.74,20.15.  Redistr  Million  subject  to  ASCC  license  01  copyrrght  seo  http://pubs.asce  01  g-'copyrrght 


Concluding  Remarks 

In  order  to  account  for  flow  variability  at  upstream  and  lateral 
boundary  conditions,  a  channel  can  be  divided  into  several  sub- 
reaches  with  representative  values  such  as  «.  T'tx).  and  K,  the  so¬ 
lutions  can  then  be  applied  to  the  individual  subreaehes.  Further, 
the  analytical  solutions  are  explicit  formulas  consisting  of  simple 
functions  with  a  number  of  controlling  parameters.  The  writers 
propose  as  future  work  to  make  an  Internet-accessible  solution 
tool,  such  as  a  small  computer  code  or  Excel  spreadsheet  to  allow 
users  to  easily  set  up  inputs  and  obtain  solutions. 
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Notation 

The  following  symbols  are  used  in  this  technical  note: 

A  =  amplitude  in  Eq.  (11); 

B  =  parameter  in  Eq.  (12); 

Cp  =  heat  capacity; 

H(t)  =  time-dependent  upstream  condition  for 
temperature; 
h  =  water  depth; 

/  =  integers,  >1; 

K  =  coefficient  of  heat  transfer; 

K  —  air/water  exchange  coefficient; 

Le  =  characteristic  envelope  length; 

LT  =  characteristic  transient  length; 
s  =  distance  downstream; 

T  =  mean  cross-sectional  temperature; 

T'bd  =  lateral  boundary  condition  for  temperature; 

Tt  —  amplitude  of  lateral  boundary  condition 
temperature  due  to  component  i; 

Tu p  =  inflow  temperature  at  an  upstream  end; 

T0  =  upstream  temperature  oscillation  amplitude, 

T\  =  lateral  boundary  oscillation  amplitude  for 
temperature; 

T\  =  downstream  gradient  of  the  amplitude  of  the 
lateral  boundary  temperature; 
f{)  =  average  upstream  inflow  temperature; 

T,  =  mean  lateral  boundary  condition  for 
temperature; 

T[  =  downstream  gradient  of  mean  lateral  boundary 
temperature; 

/  =  time, 

U  =  constant  flow  velocity; 
u  —  cross-sectional  mean  flow  velocity; 

V%y)  =  variable  flow  velocity; 

a  =  phase  of  temperature  oscillation  at  an 
upstream  end; 


p,  =  phase  of  lateral  boundary  condition  oscillation 
for  temperature  due  to  component  /; 

7,k  =  phases  in  Eqs.  (11)  and  (12); 

0,  =  phases  in  Eq.  (18); 

0,  =  phase  in  Eq.  (10); 

0J  =  phase  in  Eq.  (17); 
p  =  water  density; 

co,  =  frequency  of  lateral  boundary  temperature 
oscillation  due  to  component  /; 
co0  =  frequency  of  temperature  oscillation  at  an 
upstream  end;  and 

co,  =  frequency  of  temperature  oscillation  at  a 
lateral  boundary. 
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